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This  work  deals  with  the  assessment  of  the  thermo-kinetic  properties  of  Mg-Fe  based  materials  for 
hydrogen  storage.  Samples  are  prepared  from  MgxFe  (x:  2,  3  and  15)  elemental  powder  mixtures  via 
low  energy  ball  milling  under  hydrogen  atmosphere  at  room  temperature.  The  highest  yield  is  obtained 
with  Mgi5Fe  after  150  h  of  milling  (90wt%  of  MgH2).  The  thermodynamic  characterization  carried  out 
between  523  and  673  K  shows  that  the  obtained  Mg-Fe-H  hydride  systems  have  similar  thermodynamic 
parameters,  i.e.  enthalpy  and  entropy.  However,  in  equilibrium  conditions,  Mgi5Fe  has  higher  hydro¬ 
gen  capacity  and  small  hysteresis.  In  dynamic  conditions,  Mgi5Fe  also  shows  better  hydrogen  capacity 
(4.85  wt%  at  623  K  absorbed  in  less  than  10  min  and  after  100  absorption/desorption  cycles),  reasonably 
good  absorption/desorption  times  and  cycling  stability  in  comparison  to  the  other  studied  compositions. 
From  hydrogen  uptake  rate  measurements  performed  at  573  and  623  K,  the  rate-limiting  step  of  the 
hydrogen  uptake  reaction  is  determined  by  fitting  particle  kinetic  models.  According  to  our  results,  the 
hydrogen  uptake  is  diffusion  controlled,  and  this  mechanism  does  not  change  with  the  Mg-Fe  proportion 
and  temperature. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  recent  years  the  demand  for  an  alternative  energy  carrier  for 
vehicular  applications  has  become  more  pronounced  due  to  the 
increasing  price  of  oil-derived  fuels  and  the  pollutant  emissions 
that  their  burning  produces  (C02,  NO*).  Hydrogen  is  one  of  the 
most  attractive  alternatives  owing  to  its  abundance,  easy  synthe¬ 
sis  and  non-polluting  nature  when  used  in  fuel  cells.  Nevertheless, 
for  its  commercial  application  one  of  the  technological  constraints 
that  must  be  tackled  is  the  storage  of  this  highly  inflammable 
gas  in  an  efficient  and  safe  way.  The  relatively  advanced  method 
such  as  high-pressure  gas  and  liquefied  hydrogen  tanks  cannot 
fulfill  the  established  storage  goals  (6wt%  H  and  45kgL-1  for 
2010)  [1],  since  they  are  neither  safe  nor  density-efficient  stor¬ 
age  procedures  (hydrogen  gas:  0.99  x  1022  H  atoms  cm-3;  liquefied 
hydrogen:  4.2  x  1022  H  atoms  cm-3).  Solid-state  storage  of  hydro¬ 
gen  in  metal  hydride  form  gives  a  safety  advantage  over  the  gas  and 
liquid  storage  methods.  In  addition,  metal  hydrides  have  higher 
hydrogen  storage  density  (6.5  x  1022  H  atoms  cm-3  for  MgH2) 
[2]. 
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Magnesium  hydride  is  considered  as  a  potential  candidate  for 
solid-state  hydrogen  storage  for  on-board  applications.  MgH2  has 
the  highest  energy  density  (9  MJ  kg-1  Mg)  of  all  reversible  hydrides 
applicable  for  hydrogen  storage.  It  also  combines  a  high  hydrogen 
capacity  of  7.6  wt%,  low  cost  and  reversibility.  Furthermore,  it  is 
abundantly  available.  However,  the  practical  use  of  MgH2  is  pre¬ 
cluded  by  the  poor  rates  of  its  hydrogen  release/uptake  reactions, 
especially  at  temperatures  below  573  K.  This  fact  is  ascribed  to  the 
strong  ionic  bond  between  Mg  and  H,  which  makes  MgH2  ther¬ 
modynamically  very  stable  within  the  practical  useful  temperature 
range  [3]. 

Many  efforts  have  been  devoted  to  improve  the  hydrogen 
storage  properties  of  magnesium,  such  as  the  production  of 
micro/nanocrystalline  powder  and  creation  of  defects  by  ball 
milling  [4-9],  synthesizing  composite  of  metal  hydrides  [10-13] 
and  adding  catalyst  [4,10,14-31].  Ball  milling  is  a  well-known 
method  for  processing  metallic  powders  which  can  reduce  the  crys¬ 
tallite  size  down  to  the  nanometer  range  and  disperse  an  added 
catalyst  effectively.  Reactive  ball  milling  (RBM)  in  hydrogen  atmo¬ 
sphere  is  a  one  step  synthesis  method  which  involves  the  in  situ 
synthesis  of  hydrides  phases  facilitated  by  the  accompanied  refine¬ 
ment  of  the  microstructure  of  the  materials  during  the  process.  In 
the  case  of  magnesium,  the  crystallite  size  reduction  and  the  proper 
dispersion  of  a  3d-transition  metal  can  significantly  improve  the  H 
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transport  properties  and  hydrogenation-dehydrogenation  kinetics 
[9,14,21,25,31]. 

One  of  the  3d-transition  elements  added  to  Mg  in  order  to 
improve  its  hydrogen  storage  properties  is  Fe.  The  interest  in  the 
Mg-Fe-H  system  is  based  on  its  better  hydrogen  storage  proper¬ 
ties  compared  to  the  Mg-FI  hydride  system.  The  former  system  has 
lower  dissociation  temperature  and  pressure  than  the  latter  sys¬ 
tem  and  reasonably  good  capacity,  5.47  wt%  [32].  Despite  the  fact 
that  Mg-Fe-FI  is  considered  a  promising  material  for  hydrogen  stor¬ 
age,  many  of  its  features,  like  hydrogen  capacity,  thermodynamic 
stability  and  thermal  conductivity,  must  be  enhanced. 

In  order  to  improve  the  hydrogen  storage  properties  of  the 
Mg-Fe-FI  system,  in  the  present  work  MgxFe  (x:  2,  3  and  15)  are 
investigated.  Such  materials  are  synthesized  via  RBM  in  hydrogen 
atmosphere  at  room  temperature.  Thermodynamics  of  the  as- 
milled  materials  is  assessed  in  the  temperature  range  of  523-673  K 
by  pressure-composition  isotherm  (PCI)  plots  (reversibility  and 
equilibrium  pressure)  and  van’t  Floff  graphs  (heat  of  hydrogen 
absorption/desorption).  The  effects  of  the  pressure  and  tempera¬ 
ture  on  the  kinetic  behavior  of  the  Mg-Fe-FI  system  are  evaluated. 
The  kinetic  behavior  is  studied  between  523  and  673  K,  compar¬ 
ing  the  hydrogen  uptake  capacity,  absorption/desorption  rates  and 
cycling  stability  of  the  synthesized  Mg-Fe  compositions.  At  573  and 
623  K,  typical  kinetic  models  are  applied  to  determine  the  rate- 
limiting  step  of  the  hydrogen  absorption  reaction.  Throughout  this 
study  we  conclude  that  the  material  with  Mgi5Fe  composition  has 
better  hydrogen  storage  properties. 

2.  Experimental 

2.1.  Preparation  of  the  Mg-Fe-H  hydride  system 

The  synthesis  of  the  Mg-Fe-H  system  was  performed  via  low 
energy  ball  milling  in  hydrogen  atmosphere,  i.e.  RBM.  Stoichiomet¬ 
ric  MgxFe  (x:  2, 3  and  15)  elemental  powder  mixtures  (purity  of  the 
starting  materials  higher  than  99%)  were  mechanically  milled  in 
a  magneto-mill  Uni-Ball-Mill  II  (Australian  Scientific  Instruments) 
device,  using  stainless  steel  balls  as  a  grinding  medium.  The  milling 
details  are  in  Table  1 . 

The  vial  was  opened  several  times  during  the  milling  processes 
in  order  to  take  samples  of  the  as-milled  powders.  Before  and  after 
subtracting  the  samples  the  vial  was  evacuated  and  purged  several 
times  with  argon  gas  (+99.9995%).  During  the  RBM  procedures  the 
vial  was  refilled  with  hydrogen  each  10  h  to  maintain  hydrogen  gas 
pressure  of  0.5  MPa.  All  handling  was  carried  out  in  a  globe  box  with 
argon  gas  and  oxygen  and  moisture  controlled  atmosphere,  so  as 
to  minimize  the  oxidation  of  the  samples. 

2.2.  Characterization 

The  microstructural  and  hydrogen  storage  properties  of  the 
samples  extracted  during  the  milling  processes  and  the  obtained 
hydride  systems  were  characterized  via  diverse  techniques.  X-ray 
diffraction  (XRD)  analysis  on  a  Philips  PW 1710/01  Instruments  with 
CuKa  radiation  (X  =  1.5405  A,  graphite  monochromator,  30  mA  and 
40  kV)  was  employed  to  observe  the  present  phases  in  the  sam¬ 
ples  and  to  determine  the  crystallite  size.  The  X-ray  intensity  was 
measured  over  diffraction  angle  20  from  10°  to  90°  with  a  scan¬ 
ning  rate  of  0.02°  s-1 .  Considering  the  crystallite  size  as  the  unique 
component  responsible  for  the  XRD  peak  broadening,  the  degree 
of  microstructural  refinement  of  MgH2  and  unreacted  Fe  of  the  as- 
milled  powders  was  estimated  applying  the  Scherrer  equation.  The 
microstructural  morphology  of  the  samples  was  studied  by  scan¬ 
ning  electron  microscopy  (SEM  515,  Philips  Electronic  Instruments) 


Fig.  1.  X-ray  diffraction  patterns:  Mg2Fe,  Mg3Fe  and  Mgi5Fe  at  30  h  and  after  the 
RBM  process. 

on  stick  dispersed  powders,  and  resin-mounted  and  polished 
samples.  The  thermal  behavior  was  investigated  by  differential 
scanning  calorimetry  (DSC)  using  a  TA  2910  calorimeter  at  a  heat¬ 
ing  rate  of  5Kmin-1  and  an  argon  flow  rate  of  122  ml  min-1.  The 
amount  of  absorbed  hydrogen  and  the  proportion  of  hydride  phases 
were  estimated  from  the  DSC  curves,  using  the  peak  area  and  the 
heats  of  formation  reported  in  this  work. 

A  modified  Sieverts-type  device  coupled  with  a  mass  flow  con¬ 
troller  was  used  to  perform  both  the  thermodynamic  and  kinetic 
characterization  [33].  Absorption/desorption  PCIs  were  measured 
in  the  temperature  range  of  523-623  K  in  order  to  study  the  ther¬ 
modynamics  of  the  Mg-Fe-H  systems.  The  equilibrium  plateau 
pressures  were  calculated  as  an  average  of  the  experimental  points 
in  the  plateau  region,  taking  into  account  the  error  propagation 
theory  to  assign  an  error  range  to  each  calculated  equilibrium 
pressure.  From  the  calculated  equilibrium  pressures  the  absorp¬ 
tion/desorption  van’t  Hoff  graphs  were  plotted. 

Hydrogen  absorption/desorption  rate  measurements  were  car¬ 
ried  out  between  573  and  623  K  in  the  Sieverts-type  device.  In 
order  to  compare  the  temperature  effects  on  the  hydrogen  uptake 
rates  and  capacities  of  Mgi5Fe,  Mg3Fe  and  Mg2Fe,  all  measure¬ 
ments  were  performed  at  the  same  relationship  between  the  initial 
pressure  at  which  the  absorption  process  starts  and  the  respec¬ 
tive  equilibrium  pressure  (Pi/Pequiiibrium  =  2.4).  This  ensures  that  the 
absorption  kinetics  is  not  influenced  either  by  different  pressure 
drops  at  the  tested  temperatures  or  by  the  equilibrium  pressure 
of  the  hydride  system.  Due  to  the  experimental  fluctuations,  each 
hydrogen  absorption  measurement  was  performed  three  times,  so 
that  the  hydrogen  capacity  was  calculated  as  an  average.  Measure¬ 
ments  of  the  hydrogen  desorption  rate  were  performed  through 
a  mass  flow  controller.  These  measurements  only  give  an  estima¬ 
tion  of  the  time  required  for  the  hydrogen  desorption,  since  the 
Mg-Fe-H  systems  under  study  release  more  hydrogen  than  the 
flow  mass  controller  can  intake. 

3.  Results  and  discussion 

3.1.  Preparation  of  Mg-Fe-H  hydride  systems  by  reactive  ball 
milling 

Mg-Fe-H  hydride  systems  were  synthesized  via  RBM  in  hydro¬ 
gen  atmosphere  using  as  initial  materials  MgxFe  (x:  2,  3  and  15); 
the  details  of  the  milling  processes  can  be  seen  in  Table  1.  In  Fig.  1 
are  presented  the  XRD  diffraction  patterns  at  30  h  and  150-200  h 
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Table  1 

Processing  parameters  for  the  reactive  mechanical  milling  carried  out  with  different  Mg-Fe  stoichiometric  mixtures. 
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a  Hydrogen  purity:  +99.999%. 

b  Mg  powder  with  an  agglomerate  size  distribution  between  130  and  260  |xm  and  Fe  powder  with  an  agglomerate  size  distribution  between  40  and  130  p,m  (estimated  by 
particle  size  distribution  analysis— Mastersizer  Micro.  MAF  5000). 
c  Mg  and  Fe  powders  with  an  agglomerate  size  >1000  pirn. 


Fig.  2.  Secondary  electron  micrographs  showing  the  conglomerate  size  distribution  (on  the  left)  and  backscattered  electrons  micrographs  showing  the  dispersion  of  Fe 
(bright  phase)  and  Mg  (dark  phase)  of  the  as-milled  powders  (on  the  right).  (A)  and  (B)  Mg2Fe  after  200  h  of  RBM,  (C)  and  (D)  Mg3Fe  after  150  h  of  RBM  and  (E)  and  (F)  Mgi5Fe 
after  150  h  of  RBM. 
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of  RBM.  At  30  h  of  RBM,  in  all  the  investigated  compositions  p- 
MgH2  (JCPDS  Powder  Diffraction  Data  Card  No.  12-0697),  Fe  (JCPDS 
Powder  Diffraction  Data  Card  No.  06-0696)  and  Mg  (JCPDS  Pow¬ 
der  Diffraction  Data  Card  No.  35-0821)  phases  were  identified. 
The  XRD  patterns  taken  between  30  h  and  the  final  milling  time 
are  not  shown  as  the  unique  hydride  phase  detected  was  p-MgH2. 
In  the  case  of  the  Mg2Fe  composition,  at  200  h  of  RBM  two  small 
peaks  of  the  intermetallic  hydride  Mg2FeFl6  (JCPDS  Powder  Diffrac¬ 
tion  Data  Card  No.  38-0843)  appeared,  coexisting  with  p-MgH2 
and  unreacted  Fe.  Considering  the  reported  information  about  the 
synthesis  of  Mg2FeFl6  via  low  energy  RBM  from  2Mg-Fe,  the  fact 
that  the  ternary  complex  was  formed  after  long  hours  of  milling  is 
attributed  to  the  insufficient  energy  input  due  to  the  low  energy 
milling  mode,  insufficient  nanostructurization  of  Fe  and  the  insuf¬ 
ficient  nanocrystalline  p-MgFI2  required  to  react  with  Fe  [34-37]. 
Moreover,  the  low  intensity  of  the  Mg2FeFl6  peaks  is  attributed  to 
the  high  amorphisation  degree  reached  after  200  h  of  RBM. 

After  milling,  in  the  XRD  pattern  of  Mg3Fe  (Fig.  1)  only  p-MgFI2 
and  unreacted  Fe  were  detected.  In  the  case  of  as-milled  Mgi5Fe, 
p-MgFI2  and  Fe  were  also  present  together  with  the  metastable  7- 
MgH2  (JCPDS  Powder  Diffraction  Data  Card  No.  35-1184)  and  MgO 
(JCPDS  Powder  Diffraction  Data  Card  No.  45-0946).  Metastable  7- 
MgH2  appeared  owing  to  the  structural  defects  and  the  mechanical 
deformation  that  occur  during  the  RBM  [6].  The  presence  of  MgO 
only  in  Mg15Fe  composition  can  be  attributed  to  a  larger  surface 
area  of  free  Mg  generated  by  the  milling  process,  large  proportion 
of  Mg  particles  which  are  not  coated  with  Fe  and  the  high  reactivity 
of  Mg  to  oxygen. 

An  estimation  of  the  crystallite  size  reduction  was  carried  out 
by  Scherrer  equation.  For  Mg2Fe,  Mg3Fe  and  Mg15Fe,  the  crystallite 
size  of  MgH2  (lattice  plane  (1  10),  20:  27.9°)  and  Fe  (lattice  plane 
(111),  20:  44.5°)  reached  stable  values  after  100  h  of  RBM  (XRD 
not  shown).  Taking  into  account  the  crystallite  size  of  the  initial 
materials  (Mg2Fe:  Fe  ~40  nm  and  Mg  ~50  nm;  Mg3Fe  and  Mgi5Fe: 
Fe  and  Mg  >1000  nm)  and  the  peak  broadening  observed  in  all  XRD 
patterns  of  the  samples  (Fig.  1 ),  it  is  clear  that  the  crystallite  size  has 
been  reduced  as  a  direct  consequence  of  the  RBM  process.  The  final 
crystallite  size  values  of  MgFI2  and  Fe  are  in  the  order  of  ~10nm 
and  ~30  nm  respectively. 

Fig.  2  shows  secondary  electrons  (on  the  left)  and  backscat- 
tered  electrons  (on  the  right)  images  of  the  as-milled  powders. 
The  agglomerate  size  distributions  are  between  10  and  150  p,m  for 
Mg2Fe  (Fig.  2A)  and  10  and  300  pun  for  Mg3Fe  and  Mg15Fe  com¬ 
positions  (Fig.  2C  and  E).  The  conglomerate  average  sizes  of  the 
as-milled  powders  range  between  20  and  30  p,m  and  present  a 
spherical  porous  kind  morphology.  This  evidences  that  the  con¬ 
glomerates  size  of  the  materials  has  been  reduced  in  comparison 
to  the  initial  materials  (see  Table  1).  Micrographs  taken  on  resin- 
mounted  and  polished  powders  show  the  intermixture  degree 
between  Mg  and  Fe.  For  Mg2Fe  composition,  the  iron  particles 
embedded  in  the  Mg  matrix  cannot  be  distinguished  showing  a 
good  intermixing  degree  and  particle  size  reduction  (Fig.  2B).  On 
the  other  hand,  particles  of  Fe  bigger  than  50  p,m  can  be  observed 
in  the  micrographs  of  Mg3Fe  and  Mg15Fe  compositions  (Fig.  2D 
and  F).  Flence,  the  intermixture  degree  and  particle  size  reduction 
reached  in  Mg3Fe  and  Mgi5Fe  is  not  as  effective  as  in  Mg2Fe.  Con¬ 
sidering  the  same  milling  conditions,  this  fact  can  be  attributed  to 
the  smaller  conglomerate  size  of  the  initial  materials  (see  Table  1) 
and  the  further  refinement  degree  rendered  by  larger  amounts  of 
Fe  [9]. 

The  reduction  of  the  crystallite  and  agglomerate  size  is  caused 
by  the  mechanical  action  of  the  grinding  medium,  the  embrittle¬ 
ment  of  Mg  by  hydrogen  gas  and  the  in  situ  formation  of  MgH2 
brittle  phase  [34].  While  the  crystallite  size  reduction  down  to  the 
nanometer  range  enhances  the  hydrogen  transport  properties  of 


Fig.  3.  Hydrogen  desorption  from  the  as-milled  Mg2Fe,  Mg3Fe  and  Mgi5Fe  compo¬ 
sitions. 

the  hydride  materials,  the  reduction  of  the  agglomerate  size  results 
in  an  increase  of  specific  area  for  hydrogen  absorption.  Both  effects 
improve  the  kinetic  characteristics  of  the  Mg-Fe  based  materials 
[9,10,23,25]. 

In  Fig.  3  are  shown  the  DSC  curves  of  the  synthesized  MgFI2-Fe 
compositions.  From  these  curves  the  maximum  desorption  temper¬ 
atures  are  subtracted  and  the  relative  amount  of  hydride  phases  and 
hydrogen  capacity  estimated  (with  an  error  of  ±5%).  The  amounts 
of  unreacted  Mg  and  Fe  are  calculated  taking  into  account  the  per¬ 
centages  of  hydride  phases  and  the  stoichiometric  proportions  of 
the  initial  materials,  see  Table  2.  The  initial  desorption  tempera¬ 
ture  is  about  475  K  and  it  does  not  change  with  the  quantity  of  Fe, 
from  7  to  39wt%.  Moreover,  as  proportion  of  Fe  increases  the  total 
hydrogen  capacity  decreases  from  6.9  to  4.1  wt%  H. 

The  non-symmetrical  shape  of  the  DSC  curve  of  Mg2Fe  can  be 
assigned  to  the  presence  of  Mg2FeFl6  identified  in  the  XRD  pattern, 
see  Fig.  1  (Section  3.1),  as  Mg2FeFl6  decomposes  at  lower  temper¬ 
atures  than  MgFI2  [37].  The  deconvolution  of  the  DSC  curve  of  the 
as-milled  powders  gives  an  estimation  of  the  relative  amounts  of 
MgFI2  and  Mg2FeFI6,  see  Table  2.  It  is  considered  that  the  amount 
of  the  ternary  hydride  obtained  from  DSC  analysis  does  not  con¬ 
tradict  the  low  intensity  peaks  observed  in  the  XRD  pattern  of 
Mg2Fe  (Fig.  1,  Section  3.1)  owing  to  the  high  degree  of  amorphi¬ 
sation  of  this  hydride  phase  after  long  milling.  The  DSC  curve  of 
Mg15Fe  also  presents  a  non-symmetrical  form.  This  is  due  to  the 
lower  decomposition  temperature  of  the  metastable  7-MgFI2  [6]. 

The  desorption  temperatures  of  the  milled  materials  are  sub¬ 
stantially  lower  than  that  of  the  non-milled  MgFI2,  ~720  K  [  16].  This 
reduction  is  a  consequence  of  two  factors.  First,  the  microstructural 
changes  induced  by  milling,  i.e.  particle  size  reduction,  crystallite 
size  reduction,  creation  of  grain  boundaries  and  defects  density 
[34,38].  Second,  the  catalytic  effect  produced  by  Fe  [34,36,38], 
since  the  desorption  temperature  of  milled  MgFI2  (milled  for  100  h, 

Table  2 

Estimation  of  the  relative  amount  of  phases  (±5%  error)  of  the  as-milled  composi¬ 
tions  and  maximum  desorption  temperature  from  DSC  analysis. 
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Table  3 

Equilibrium  pressures  and  relative  amount  of  hydrogen  stored  in  the  Mg-Fe-H  hydride  system  as  a  function  of  the  temperature,  obtained  from  the  absorption/desorption 
PCI  curves. 


Sample 

Temperature  [I<] 

Absorption 

MgH2-Mg2FeH6 

Desorption 

MgH2 

Mg2FeH6 

Peq  [kPa] 

Hydrogen  content  [wt%] 

Feq  [kPa] 

Hydrogen  content  [wt%] 

Peq  [kPa] 

Hydrogen  content  [wt%] 

Mg2Fe 

673 

1925  ±  100 

3.90 

1565  ±  30 

0.93 

900  ±  60 

2.97 

648 

1160  ±  70 

2.80 

1000  ±  40 

0.90 

470  ±  130 

1.90 

623 

730  ±  20 

2.60 

590  ±  20 

1.95 

260  ±  65 

0.65 

573 

250  ±  30 

2.00 

163  ±  11 

2.00 

- 

- 

523 

60  ±4 

1.90 

50  ±  1 

1.58 

- 

- 

Mg3Fe 

673 

2360  ±  160 

4.20 

1955  ±  45 

2.80 

915  ±  40 

0.90 

648 

1440  ±  80 

4.00 

1210  ±  40 

3.42 

510  ±  70 

0.53 

623 

865  ±  65 

3.90 

710  ±  30 

3.80 

- 

- 

573 

290  ±  60 

3.70 

200  ±  10 

3.70 

- 

- 

523 

75  ±  30 

2.80 

40  ±4 

2.80 

- 

- 

Mgi5Fe 

673 

1955  ±  43 

6.15 

1700  ±  57 

5.95 

900  ±  40 

0.20 

663 

1600  ±  35 

6.05 

1380  ±  40 

5.95 

650  ±  60 

0.10 

648 

1245  ±  80 

6.00 

1010  ±  37 

6.00 

510  ±  90 

0.05 

633 

980  ±  80 

6.00 

750  ±  24 

6.00 

- 

- 

623 

800  ±  88 

6.00 

580  ±  30 

6.00 

- 

- 

573 

190  ±  8 

5.50 

170  ±  7 

5.50 

- 

- 

548 

150  ±  48 

5.30 

88  ±6 

5.30 

- 

- 

523 

70  ±  24 

3.90 

47  ±  3 

3.90 

- 

- 

crystallite  size  of  12  nm)  is  655 1<  [39,40].  As  it  can  be  noticed, 
the  decomposition  temperature  lowers  as  the  Fe  proportion  is 
increased.  However,  a  small  quantity  of  Fe  is  enough  to  reduce  the 
decomposition  temperature  and  to  reach  better  hydrogen  capacity. 

The  reactive  ball  milling  method  has  been  used  to  synthesized 
hydride  materials  from  MgxFe  (x:  2,  3  and  15).  Via  this  one  step 
synthesis  procedure  nanostructurized  materials  composed  by  dif¬ 
ferent  amounts  of  hydride  phases  have  been  obtained.  Among  the 
investigated  compositions,  Mgi5Fe  has  the  higher  yield  of  hydride 
phase  and  hydrogen  capacity,  see  Table  2. 

3.2.  Thermodynamic  characterization 

The  equilibrium  properties  of  the  synthesized  Mg-Fe-H  systems 
were  examined  between  523  and  673  K.  Pressure-composition 
isotherms  were  measured  to  determine  the  absorption/desorption 
equilibrium  pressures  and  hydrogen  capacities  of  the  Mg2Fe,  Mg3Fe 
and  Mgi5Fe  compositions.  The  equilibrium  pressures  and  hydrogen 
reversible  capacities  were  estimated  from  the  plateau  widths.  As  it 


can  be  noticed  in  Table  3  as  well  as  in  Fig.  4,  hydrogen  capacities 
and  equilibrium  pressures  decrease  as  the  temperature  falls.  This 
is  due  to  the  atom  diffusion  mechanisms  which  are  hampered  by 
the  lower  temperatures  [41  ]. 

The  PCI  of  the  extreme  compositions,  viz.  Mg2Fe  and  Mgi5Fe, 
can  be  observed  in  Fig.  4.  The  hydrogen  absorption/desorption  sto¬ 
ichiometric  reactions  can  be  seen  in  Ref.  [42].  One  plateau  is  only 
observed  during  the  absorption  processes  because  the  equilibrium 
pressures  of  Mg2FeH6  and  MgH2  are  similar  [41].  On  the  other 
hand  the  equilibrium  pressures  of  MgH2  and  Mg2FeH6  during  the 
release  of  hydrogen  are  different,  thus  hydrogen  desorption  pro¬ 
cesses  show  two  plateaus.  In  the  case  of  Mg2Fe,  two  plateaus  can  be 
noticed  during  the  desorption  process  at  temperatures  above  623  K. 
This  is  a  result  of  the  Mg2FeH6  formation  via  thermal  enhanced 
atom  diffusion  mechanism  during  the  hydrogen  absorption.  For 
Mg3Fe  and  Mgi5Fe,  two  desorption  plateaus  are  observed  at  tem¬ 
peratures  above  648  K,  but  they  are  not  seen  at  623  K.  In  this  case, 
the  amount  of  available  Fe  is  not  in  the  right  stoichiometric  pro¬ 
portion  for  the  Mg2FeH6  formation,  the  Fe  conglomerates  are  in 


Hydrogen  Capacity  [wt%]  Hydrogen  Capacity  [wt%] 


Fig.  4.  Pressure-composition  isotherms  for  (A)  Mg2Fe  and  (B)  Mgi5Fe  compositions  at  573,  623  and  673 1<. 
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Fig.  5.  (A)  Absorption  and  (B)  desorption  van’t  Hoff  plots  of  the  (—  fitting,  A  experimental)  Mg2Fe  and  ( - fitting,  •  experimental)  Mgi5Fe  compositions. 


the  order  of  50  jxm  and  are  not  well  dispersed  as  in  Mg2Fe  compo¬ 
sition  (see  Fig.  2— backscattered  electrons).  Therefore,  below  648  K 
the  thermal  enhanced  atom  diffusion  mechanism  to  form  Mg2FeH6 
during  the  hydrogen  absorption  process  might  be  hampered  by  the 
poor  Fe-Mg  contacting.  The  presence  of  one  absorption/desorption 
plateau  below  623  K  (Mg2Fe)  and  648  K  (Mg3Fe  and  Mgi5Fe)  is 
attributed  to  MgFI2  which  is  the  unique  hydride  phase  formed  dur¬ 
ing  the  hydrogen  absorption  process  [40]. 

Another  important  characteristic  shown  in  Fig.  4  is  the  phe¬ 
nomenon  of  hysteresis  in  the  absorption/desorption  cycles.  As  the 
temperature  increases  it  becomes  more  noticeable.  It  can  be  clearly 
observed  in  the  PCIs  of  the  Mg2Fe  composition,  even  at  573  K.  On 
the  contrast,  it  is  not  too  pronounced  in  the  PCI  cycles  of  the  Mgi5Fe 
composition.  The  phenomenon  of  hysteresis  is  ascribed  to  elastic 
and  plastic  strain  energies  related  to  the  volume  changes  during 
the  formation  of  the  hydride  phases  [43].  Moreover,  at  tempera¬ 
tures  above  623  K,  at  which  Mg2  FeHg  is  formed  during  the  hydrogen 
absorption  process,  the  hysteresis  arises  manly  from  the  presence 
of  Mg2FeH6  [44,45]. 


The  thermodynamic  parameters  (absorption/desorption 
enthalpy  and  entropy)  of  Mg2Fe,  Mg3Fe  and  Mg15Fe  composi¬ 
tions  were  calculated  from  the  van’t  Floff  graphs.  In  Fig.  5  are 
shown  the  absorption/desorption  van’t  Floff  plots  of  the  extreme 
compositions,  i.e.  Mg2Fe  and  Mgi5Fe,  in  the  temperature  range 
of  523-673  K.  For  the  sake  of  clarity  and  due  to  the  similarity  of 
the  calculated  thermodynamic  parameters  of  Mg3Fe  (Table  4),  the 
van’t  Floff  plots  of  this  intermediate  composition  are  not  included 
in  Fig.  5. 

In  Table  4  the  values  of  enthalpy  and  entropy  are  shown  and 
compared  with  those  found  in  the  literature  and  our  previous  works 
[32,40-42,46-48].  The  calculated  parameters  show  a  good  correla¬ 
tion  with  those  found  in  the  literature  for  the  hydride  mixture.  The 
differences  between  the  values  of  calculated  parameters  and  those 
reported  in  the  literature  might  be  related  to  the  measurement 
procedure  and  the  temperature  ranges  considered. 

The  assessed  compositions,  i.e.  Mg2Fe,  Mg3Fe  and  Mgi5Fe, 
present  thermodynamic  parameters  that  lay  in  the  same  range.  This 
shows  that  the  thermodynamic  of  the  Mg-Fe-H  system  is  not  mod- 


Table  4 

Formation  and  decomposition  enthalpies  and  entropies  obtained  from  the  van’t  Hoff  plot  of  the  Mg-Fe-H  hydride  system  synthesized  from  different  stoichiometric 
compositions. 


Absorption 

MgH2  and  Mg2FeH6 

AH  [kj  mol-1  H2] 

AS  [J  moF1  H2  K"1  ] 

Mg2Fe  (523-673 1<) 

-68  ±  1 

-124  ±  2 

Mg3Fe  (523-673  K) 

-68  ±  5 

-127  ±  7 

Mg15Fe  (523-673  K) 

-65  ±2 

-122  ±  4 

Ref.  [41]  (548-723 1<) 

-66  ±2 

-124  ±  3 

Ref.  [48]  (698-723  K) 

-55  ±  3 

- 

Desorption 

MgH2 

Mg2FeH6 

AH  [kj  mol-1  H2] 

AS  [J  mol"1  H2  K1  ] 

AH  [kj  mol-1  H2] 

AS  [J  mol1  H2  K"1  ] 

Mg2Fe  (523-673 1<) 

-68  ±2 

—124  ±3 

-87  ±5 

-147  ±7 

Mg3Fe  (523-673  K) 

-76  ±1 

-138  ±  2 

-85 

-144 

Mg15Fe  (523-673  K) 

-71  ±1 

-128±2 

-96  ±24 

-161  ±36 

Ref.  [41  ]  (548-648  K)1,  (573-648  K)2 

-67±21 

—  123  A31 

-80  ±72 

-137  ±132 

Ref.  [42]  (573-673  K) 

-79  ±3 

-140  ±7 

-87  ±3 

-147  ±15 

Ref.  [46]  (587-849  K) 

-74 

135 

- 

- 

Ref.  [47]  (648-723  K) 

- 

- 

-98  ±3 

- 

Ref.  [48]  (623-698  K) 

—77  ±4 

—138  ±  3 

-86  ±6 

-147  ±9 

Ref.  [32]  (623-798 1<) 

- 

- 

-77 

- 
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Table  5 

Hydrogen  absorbed  in  500  s  at  523,  573,  623  and  673  K  after  a  number  of  hydrogen  absorption/desorption  cycles,  between  parentheses. 


Sample 

Theoretical  capacity  H  [wt%] 

H  [wt%]  absorbed  in  500  s  at  different  temperatures 

523  K 

573  K 

623  K 

673  K 

Mg2Fe 

3.8 

0.95  (80) 

1.35(100) 

1.70(110) 

2.50(120) 

Mg3Fe 

4.7 

1.40(10) 

2.4(20) 

3.45  (35) 

3.50(50) 

Mgi5Fe 

7.1 

1.50(50) 

3.00  (60) 

4.85(110) 

4.90(100) 

Table  6 

Integral  forms  of  various  gas-solid  kinetics  models  used  for  fitting  experimental  sorption  data  [49,50]. 


Model 

Description 

Integral  form  g(a)  =  k  t 

Nucleation  models 

Johnson-Mehl-Avrami  (JAM),  n  =  2 
Johnson-Mehl-Avrami  (JAM),  n  =  3 

Two-dimensional  growth  of  existing  nuclei  with  constant  interface  velocity. 
Three-dimensional  growth  of  existing  nuclei  with  constant  interface  velocity. 

[ln(l-a)]1/n 

Geometrical  contracting  models 

Contracting  area  (CA) 

Contracting  volume  (CV) 

Two-dimensional  growth  with  phase  boundary  controlled  reaction. 

Three-dimensional  growth  with  phase  boundary  controlled  reaction. 

1  -  (1  -  a)1/2 

1  — (1  -a)1/3 

Diffusion  models 

ID  diffusion 

2D  diffusion 

3D  diffusion  (Ginstling-Brounshtein) 

Surface  controlled  (chemisorption). 

Two-dimensional  growth  diffusion  controlled  with  decreasing  interface  velocity. 
Three-dimensional  growth  diffusion  controlled  with  decreasing  interface  velocity. 

O'2 

[(1  -a)ln(l  -a)\  +  a 

1  -  (2/3 )a  -  (1  -  a)2/3 

ified  by  the  stoichiometric  proportions  of  Mg  and  Fe.  It  is  also  worth 
noting  that  Mgi5Fe  has  the  advantages  of  higher  capacity  than 
the  other  compositions  and  small  hysteresis  between  the  hydro¬ 
gen  uptake  and  release.  These  two  features  are  quite  important  for 
potential  practical  applications. 

3.3.  Kinetics:  hydrogen  realize  I  uptake  behavior 

The  kinetic  behavior  was  investigated  between  523  and  673  K 
using  an  established  relationship  of  P/Pequiiibrium  equal  to  2.4  in 


a  Sieverts-type  device.  Measurements  of  the  hydrogen  absorp¬ 
tion  capacity  reached  in  500  s  and  after  different  numbers  of 
absorption/desorption  cycles  are  shown  in  Table  5.  The  reported 
capacities  correspond  to  MgH2  in  presence  of  free  Fe.  At  an  inter¬ 
mediate  temperature,  623  K,  the  hydride  material  composed  by 
Mg15Fe  shows  a  capacity  of  4.85  wt%  H  absorbed  in  less  than  10  min 
and  after  100  absorption/desorption  cycles.  These  hydrogen  stor¬ 
age  uptake  characteristics  are  better  in  comparison  to  the  other 
compositions  and  are  also  interesting  conditions  for  a  potential 
application. 


Time  [s]  Time  [s] 


Fig.  6.  Hydrogen  absorption/desorption  rate  for  the  different  Mg-Fe  compositions:  (A)  absorption  kinetics  and  (B)  desorption  kinetics  at  623  K;  (C)  absorption  kinetics  and 
(D)  desorption  kinetics  at  573  K. 
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Fig.  7.  Curves  resulted  from  the  application  of  the  integral  forms  of  the  particle  kinetics  models  to  the  hydrogen  absorption  rate  measurements  at  573  and  623  K  for  the  different 
Mg-Fe  compositions  (Fig.  6).  —  Linear  fits;  the  absorbed  hydrogen  fraction,  H.F.,  followed  by  the  range  of  time  required  to  reach  the  specified  H.F.  (between  parentheses)  and 
the  correlation  coefficient,  R,  are  indicated. 


A  comparison  between  the  hydrogen  absorption  and  desorp¬ 
tion  rates  at  523  and  623 1<  is  shown  in  Fig.  6.  At  623 1<  (Fig.  6A),  it  is 
clear  that  Mgi5Fe  is  faster  than  the  other  compositions.  However, 
at  573 1<  (Fig.  6C),  the  rate  of  hydrogen  uptake  of  Mg3Fe  is  slightly 
faster  than  Mgi5Fe.  The  hydrogen  release  behavior  of  the  hydride 
materials  is  shown  in  Fig.  6B  and  D.  Due  to  the  experimental  limita¬ 
tions  (see  Section  2.2),  these  measurements  are  an  estimation  of  the 
times  required  for  the  hydrogen  desorption.  At  573  and  623  K,  des¬ 
orption  times  range  between  400  and  600  s.  In  the  case  of  Mgi5Fe, 
which  is  the  composition  with  higher  capacity,  at  573  and  623  K  the 
hydrogen  release  takes  400  and  500  s  respectively. 

In  order  to  obtain  information  about  the  reaction  mechanism 
during  the  hydrogen  uptake,  an  analysis  of  the  curves  shown  in 
Fig.  6A  and  C  applying  gas-solid  kinetic  models  was  performed. 
The  general  expression  to  study  the  rate  of  gas-solid  reactions  is  as 
follows 

^ =K(T)xF(P)xC(a )  (1) 

where  the  overall  reaction  rate  is  function  of  the  temperature 
(T),  hydrogen  gas  pressure  (P)  and  the  absorbed  hydrogen  frac¬ 
tion  (a— defined  in  this  work  as  fraction  of  the  total  hydrogen 
wt%  reached  in  500  s).  The  temperature  term  is  a  function  of  the 
Arrhenius  equation  I<(T)  =  A  x  e~(EalRT)t  and  the  F(P)  term  is  taken  as 


^/^equilibrium  equal  to  2.4.  Considering  that  JC(T)  does  not  change  at 
given  temperature  and  F(P)  is  constant,  Eq.  (1 )  can  be  expressed  as 

*«>-/”< SSj-**'  121 

where  g(o t)  is  the  integral  form  of  the  models,  which  are  shown  in 
Table  6. 

Although  metal-hydrogen  reactions  involve  a  complex 
sequence  of  steps,  in  the  above  context  a  simplified  scheme 
with  three  steps  can  be  proposed.  In  this  sense,  hydriding  reaction 
involves  a  first  stage  starting  at  the  gas-surface  region,  following 
by  the  nucleation/growth  of  the  hydride  phase  and  ending  with 
the  diffusion  through  the  transformed  phase.  When  absorption 
rate  is  governed  by  one  of  these  steps,  the  kinetic  rate-limiting  step 
can  be  identified  by  fitting  the  integral  equations  of  the  particle 
reaction  models  (Table  6). 

The  reaction  of  Mg  and  Mg-based  materials  with  hydrogen  is 
influenced  by  three  dominant  factors:  crystallite  size,  particle  size 
and  catalyst  [50].  From  the  application  of  the  integral  forms  of  var¬ 
ious  gas-solid  kinetics  models  (Table  1 )  to  the  absorbed  hydrogen 
fraction,  a ,  reached  in  500  s  (Fig.  7),  it  is  evident  that  the  mechanism 
is  neither  surface  controlled  nor  nucleation/growth  controlled.  This 
behavior  is  attributed  to  the  catalytic  effect  of  Fe  on  the  hydro- 
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gen  dissociation  and  the  refinement  degree  of  the  microstructure, 
respectively.  As  it  can  be  observed,  for  all  the  practised  composi¬ 
tions  the  best  fit  is  obtained  with  the  3D  diffusion  model,  which 
indicates  that  the  hydrogen  mechanism  is  controlled  by  the  diffu¬ 
sion  of  hydrogen  into  the  MgH2  bulk.  Our  results  are  in  agreement 
with  the  previously  reported  for  hydrogen  absorption  in  milled 
MgH2  catalyzed  with  Nb203  [50].  Then,  the  microstructural  refine¬ 
ment  and  the  dispersion  of  Fe  imparted  by  RBM  have  enhanced  both 
hydrogen  dissociation  and  nucleation/growth  rate,  without  effect 
on  the  hydrogen  diffusion  coefficient  through  the  hydride  phase. 

4.  Conclusion 

Mg-Fe  based  materials  for  hydrogen  storage  were  produced  via 
RBM  using  mixtures  of  MgxFe  (x:  2,  3  and  15)  elemental  powders. 
Via  this  one  step  method  of  synthesis,  a  high  degree  of  microstruc¬ 
tural  refinement  and  catalyst  dispersion  (Fe),  which  enhanced  the 
hydrogen  storage  properties,  was  reached.  The  higher  yield  of 
hydride  phase  (90  wt%  of  MgFI2 )  and  hydrogen  capacity  (6.9  wt%  FI) 
was  obtained  from  the  Mgi5Fe  composition.  Studies  under  equilib¬ 
rium  conditions  showed  that  the  thermodynamic  parameters  for 
different  Mg-Fe  compositions  are  similar  and  are  in  good  agree¬ 
ment  with  those  reported  in  the  literature.  An  assessment  of  the 
hydrogen  uptake  reaction  mechanism  by  fitting  kinetic  models 
showed  that  the  absorption  reaction  is  diffusion  controlled.  This 
mechanism  does  not  change  with  the  Mg-Fe  proportion  at  the 
measured  temperatures. 

From  the  studied  compositions,  Mg15Fe  showed  the  highest 
hydrogen  capacity  (4.85  wt%  H  after  100  absorption/desorption 
cycles  at  623  K),  reasonable  absorption/desorption  times  (absorbs 
90%  of  its  final  capacity  in  about  200  s  and  the  hydrogen  release 
takes  600  s)  and  narrow  hysteresis  cycle.  These  are  promising  char¬ 
acteristics  from  the  application  point  of  view. 
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